import sympy

x=sympy.Symbol('x')
f=sympy.pi*x**2/2
res=sympy.integrate(f,(x,0,4))
sympy.pprint(res)